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Abstract 

A potential model for four interacting quarks is constructed in SU(2) from six basis states - the three partitions into 
quark pairs, where the gluon field is either in its ground state or first excited state. With four independent 
parameters to describe the interactions connecting these basis states, it is possible to fit 100 pieces of data - the 
ground and first excited states of configurations from six different four-quark geometries calculated on a 16^ x 32 

lattice. 
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I. INTRODUCTION 



Over the years there has been much progress in lattice QCD. However, this has been restricted to systems with 
only a few quarks - upto three in most cases. On the other hand, in particle and nuclear physics there is also 
considerable interest in few- and miilti-hadron systems beginning with the possibility of bound KK states. It seems 
unlikely that direct results for any but the simplest of these hadron systems could be obtained in the near future from 
lattice QCD simulations. Our approach is thus to construct a model that reproduces lattice results for the simplest 
multi-quark systems in a way that can be relatively easily extended to more complex cases. The lattice data has 
been obtained earlier by the Helsinki group [l](and references therein) as energies of four static quarks in various 
geometries, such as quarks at the corners of squares, rectangles, tetrahedra and some other less symmetric geometries. 
This 'experimental' data is then to be explained by a model. The choice of geometries is hopefully general enough for 
the model to reproduce also the energies for geometries intermediate to those actually considered. 

It should be pointed out that the 'experimental' data is not for full QCD. For numerical reasons, the SU(2) gauge 

group is used instead of SU(3), the use of the latter requiring an order of magnitude increase in computer resources. 
However, the indications arc that, for the present type of discussion, SU(2) suffices. Another approximation is taking 
the quarks to have an infinite mass. This means that our constituent quarks are static and that there are no sea 
quarks (the so-called quenched approximation). At present, this limit is now being partially removed by applying the 
techniques of Ref. [2] to a system of two light and two heavy (static) quarks. The sea quarks will appear through 
using gauge configurations generated by the UKQCD collaboration for full QCD. However, in this paper, we only 
attempt to understand systems with four static SU(2) quarks in the quenched approximation. 

Any model that can be extended to multi-quark systems must presumably depend explicitly only on the quark 
degrees of freedom, with the gluon degrees of freedom entering only implicitly. This is the same philosophy that has 
been successful for interacting multi-nucleon sytems. In the latter, the meson fields generating the interactions are 
'summarized' as internucleon potentials, which take the form of two-nuclcon potentials and, to a lesser extent, three- 
and four-nucleon potentials. However, an important difference between multi-quark and multi-nucleon systems is the 
strongly self-interacting nature of the gluon fields mediating the interaction in the former case. It is, therefore, not at 
all clear that any effective model defined in terms of only quark degrees of freedom will be successful. It is the purpose 
of this article to see to what extent a model can be developed for the four-quark case. If this attempt fails, then 
there is no point in expecting it to work in even more complex quark problems; a successful model for four quarks is 
necessary but not sufficient before considering any extension to even more quarks. 

It should be added that, in spite of the problems outlined above, many groups still consider that multi-quark systems 
can be treated using simply the basic two-quark potential e.g. [3,4]. In the opinion of the authors - and backed up 
by their earlier works - there seems to be no justification for such an approach. 

In Section 2 three versions of the model, starting with the simplest, are introduced and in Section 3 results are given. 



Since the four quarks in the lattice calculation are static, the corresponding model should not contain any kinetic 

energy. Also, because we only consider SU(2), there is no distinction between the group properties of quarks(g) and 
antiquarks(g). Four such quarks can then be partitioned as pairs in three different ways 



where each (qigj) is a colour singlet. However, these three basis states are not orthogonal to each other. Also, 
remembering the fact that the quarks are indeed fermions gives, in the weak coupling limit, the condition in the 
appendix of Ref. [5] 



II. THE MODEL WITH 2, 3 OR 6 BASIS STATES 



^ = {qiq3){q2q4), B = {qiqi){q2qi) and C = (91^2) («3g4), 



(2.1) 



\A + B + C>=Q. 



(2.2) 



Since < A\A >=< B\B >=< C\C >= 1, we get - in this limit - the equalities 
< A\B >=< B\C >=< A\C >= -1/2. 
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A. 2 basis states 



N=( ' /o 7'/' )• (2-3) 



Restricting the basis to only the two states A and B leads to the normalisation matrix 

-1/2 1 

In addition to this there is the potential matrix 

Y^d^A Vab\ ^2 4) 

V Vba vbb J ^ ' 

In the weak coupling limit - when all quarks are close to each other - we expect the potential matrix to simply become 
that form appropriate for one-gluon exchange i.e. 

"l^y = -^X!^ ■'^j^y' (2.5) 

where Vij = —e/rij. Then Vaa = vi3 + V2i etc. and 

Vab = Vba = ~ {vi3 + V2i + vu + ^23 - V12 - ^34) • (2.6) 

The object is then to compare the eigenvalues of 

[V - £;(4)N] = (2.7) 

with the lattice results - the success or failure of the model being to what extent the two agree. Here E(A) is the total 
four-quark energy. However, in practice, it is more convenient to deal with the four-quark binding energy E defined 
as, E = E{4) — VaA: assuming that the A partition has the lowest energy. A perturbative calculation to 0(o?) has 
reproduced this two-state model [6] . 

If all three basis states arc included, the model has a problem since the matrix in Eq. 2.7 is singular for the obvious 
reason that |j4 + iJ + C>=0. In some of our earlier work this was interpreted to mean that it was unnecessary 
to include all three states and so the symmetry was broken by keeping the two states with the lowest energy, let us 
say , A and B. A similar thing also occurred in the lattice simulations. There it was found that the energy of the 
lowest state was always the same in both a 2x2 and 3x3 description, providing A ox B had the lowest energy. In 
addition the energy of the second state was, in most cases, more or less the same - the largest difference occurring 
with the tetrahedral geometry. Another modification to the basic model in Eqs. 2.5-2.7 is to replace the original 
one-gluon-exchange potential by the full two-quark potential 

% = -e/ry + hsTij + c. (2.8) 



This is now a move away from the extreme weak coupling limit. In principle, everything is now known, since the 
parameters e,bs,c in Vij can be determined from the same lattice calculation as the four-quark energies. However, 
this model fails badly, since it still contains the van der Waals long range forces. The main feature of this failure is 
that the lowest eigenvalue (Eq) is too low and the other (Ei) too high. This suggests that the off-diagonal potential 
Vab is too large. However, as we go away from the weak coupling limit, it is expected that < A\B > will decrease 
until eventually, in the extreme strong coupling limit, it would vanish. This, therefore, suggests that for some realistic 
intermediate coupling the off-diagonal normalisation matrix element Nab should also be reduced. This effect we 
simulate by introducing a factor /, which decreases as the distance between the quarks increases i.e. 

< A\B >= -f/2. (2.9) 

For internal consistency, this same factor must also be introduced into Vab i-e. 

Vab ^ fVAB (2.10) 
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otherwise the eigenvalues would depend on the self-energy term c in Eq. 2.8. In Ref. [1] this model had the good 
feature that, when fitting the data (i^oj Ei) for a given square, only a single / was necessary to fit both energies. Of 
course, / was dependent on the size of the square, but a reasonable parametrization was 

f{Ia) = eM-b.kfS), (2.11) 



where bg was the string energy of Eq. 2.8, S the area of the square and kf w 0.5. The form of this parametrization 
was motivated by strong coupling ideas [7-9]. The original hope was that, with kf determined from the squares, the 
model should automatically also fit other geometries with 5* being the appropriate area contained by the four quarks. 
When the four quarks lie in a plane, the definition of S is clear. However, in non-planar cases the situation is more 
complicated and so here the area is simply taken to be the average of the sum of the four triangular areas defined 
by the positions of the four quarks i.e. the faces of the tetrahedon. For example, in the notation of Eq. 2.1, the 
appropriate area S{AB) for / is 

S{AB) = 0.5[S'(431) + 5(432) + S'(123) + 5(124)], (2.12) 



where S{ijk) is the area of the triangle with corners at i,j and k. For planar geometries this simply reduces to the 
expected area. But for non-planar cases this is only an approximation to S{AB) - a more correct area being one that 
is not necessarily a combination of planar areas but of curved surfaces with minimum areas. These possibilities are 
discussed in Ref. [10]. It would be feasible to incorporate this refinement here, since only a few (« 50) such areas 
are needed. But for a general situation, in which the positions of the quarks are integrated over, it would become 
impractical, since the expression for the minimum area itself involves a double integration. In contrast, the area used 
here is an algebraic expression and is, therefore, more readily evaluated for any geometry. The above model will be 
referred to as Version la. 

This model has only one free parameter kf in /. Another possibility with additional parameters /o, kp is 

f{Ib) = foexp{-bskfS + y/b'skpP) (Version lb), (2.13) 



where P is the perimeter bounding S. However, as shown in Ref. [11], this reduces in the continuum limit to the same 

as Version la - the differences at /3 = 2.4 being mainly due to lattice artefacts. 

Unfortunately, both of these models have the feature that, for regular tetrahedra, they are unable to reproduce a 
degenerate ground state with a non-zero energy, since the two eigenvalues are 

^0 = - YTJJ2 1^^^ ~ ^^^^ ''''^ " t^^^ ~ ^^^^ 



and for regular tetrahedra Vcc = Vaa- 



B. 3 basis states 



In the previous subsection only two of the possible three basis states A, B, C were used, the original reason for this 
being the condition in Eq. 2.2. However, once the factor / is introduced as in Eqs. 2.9, 2.10 this condition no longer 
holds, so that all three states can be incorporated. This apparently leads to the need for two more factors /', /" 
defined by 

<A\C>=-f'/2 and < B\C>=-f"/2. (2.15) 



However, with the parametrization of / as in Eq. 2.11 or 2.13 and the definition of S as in Eq. 2.12. it is seen that 
/' = /" = /, since S is simply proportional to the area of the faces of the tetrahedron defined by the four quark 
positions and is independent of the state combination used. Therefore, a 3 x 3 model has the form where the N and 
V matrices are: 

/ 1 -//2 -//2 \ 
N(/) = -f/2 1 -f/2 and (2.16) 
V -f/2 -f/2 1 J 
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(VAA /Vab /Vac \ 
JVba vbb JVbc 
fVcA fVcB VCC J 



(2.17) 



Unfortunately, this niodol(II) also has some unpleasant features: 

i) Again for regular tetrahedra all three eigenvalues are zero. 

ii) For a linear geometry, since the 'appropriate' area as defined by Eq. 2.12 vanishes, we get / = 1 i.e. we are back 
to the weak coupling limit and a singular matrix. 

iii) For squares the model gives Ei = —Eq, whereas the predictions of Model I in Eqs. 2.14 seem to be nearer the 
data. 

The regular tetrahedron result is clear. There is only one energy scale in the model, since all the Vij are the same. 
Therefore, there can not be any excitations. The positive feature of this model is that all three states are treated 
equally - a point that seems to be necessary for regular tetrahedra. 



The above models both have trouble in describing regular tetrahedra. This might be considered a minor point to 
worry about, since such a configuration is very special. However, our philosophy is that, if any configuration cannot 
be fitted, then the model fails, since then there is no reason to expect configurations not checked explicitly to be fitted. 
An interesting feature of the regular tetrahedron data is that the lowest state becomes more bound as the tetrahedron 
increases in size with the magnitude of Eq increasing from -0.0202(8) to -0.028(3) as the S cube containing the 
tetrahedron increases from ri = 2 to ri = 4. This is opposite to what happens with squares, where the magnitude of 
Eq decreases from —0.0572(4) to —0.047(3) as d increases from 2 to 5. This indicates that there is coupling to some 
higlic;r state(s) that becomes more effective as the size increases and suggests that these higher states contain gluon 
excitation with respect to the A,B,C configurations. Model III, introduced in Ref. [11], therefore, extends Model II 
by adding three more states A*,B*,C* , where in analogy with Eq. 2.1 



Here (giqa)^;^ denotes a state where the gluon field is excited to the lowest state with the symmetry of the E^ 
representation of the lattice symmetry group D4h. Because it is an odd parity excitation, A*,B*,C* must contain 
two such states in order to have the same parity as A, B, C. The excitation energy of an Eu state over its ground 
statc(^ig) counterpart is » tt/ R for two quarks a distance R apart. As R increases this excitation energy decreases 
making the effect of the A*,B*,C* states more important, leading to the effect mentioned above. Here we have 
assumed that these states arise from a combination of excited states with £„. However, it is possible that they involve 
other excitations, e.g. 



where the A[g state is a gluonic excitation with the same quantum numbers as the ground state (Aig). For this 
case the following formalism would be essentially the same. Another possibility, which is not considered here, is 
that the relevant excitations are fiux configurations where all four quarks, instead of two, are involved in forming a 
colour singlet. In the strong coupling approximation such states would reduce to two-body singlets due to Casimir 
scaling of the string tensions - the string tension for a higher representation would be more than double the value 
of the fundamental string tension, thus preventing junctions of two strings in the fundamental and one in the higher 
representation. This would happen both in SU(2) and SU(3), the only exception being the unexcited C state in SU(3), 
which would involve an antitriplet string. 

There are now several new matrix elements that need to be discussed for N(f) and V(f): 

a) With the addition of the A*,B*,C* states and the antisymmetry condition A* +B* +C* =0 analogous to Eq. 2.1, 
there are now two more functions f"''^ defined as 



C. 6 basis states 



A* = {qiq3)E^{q2q4)E„ etc. 



(2.18) 



A* = {qiq3)A{^{q2q4) etc. 



(2.19) 



< A*\B* >=< A*\C* >=< B*\C* >= -/72 and 
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< A*\B >=< A*\C >= .. etc. .. = -/"/2. (2.20) 

Here it is assumed that /"'^ are both dependent on S as defined in Eq. 2.12. Since involves only the excited states, 
it is reasonable to expect it has a form similar to / in Eq. 2.11 i.e. 

r = eM-bskcS). (2.21) 



b) By orthogonality < A\A* >=< B\B* >=< C\C* >= 

c) In the weak coupling limit, from the A* + B* + C* = condition, we expect < A\B* >=< B\C* >= = at 

small distances. To take this into account we parametrize /" as 

r = ifi + bsf^S) exp{-bskaS). (2.22) 

A partial justification of this is given in the Appendix and in Fig 6 of Rcf. [11]. In the following fit all three parameters 
/f , , ka are varied. However, it is found that /f is always consistent with zero - as expected from the above condition 
< A\B* >= .... = 0. Therefore, in most of this work /f is fixed at zero. 

d) For the potential matrix V(f) the diagonal matrix elements, after the lowest energy Vdd amongst the basis states 
is removed, are 

<A*\V- Vdd\A* >= ^;*(13) + v*{24) - Vdd, etc., 

where v*{ij) is the potential of the Eu state - a quantity also measured on the lattice along with the four-quark 

energies and 

Vdd = mm[VAA = «(13) + i;(24), Vbb = v{U) + w(23), Vcc = v{12) + i;(34)]. 
However, to allow more freedom in the following fits, we introduce a parameter bo, where 

<A*\V-VDD\A*>=VAA-VDD+boVXA , <B*\V-Vdd\B* >=VBB-VDD+boV^B, 

<C*\V- Vdd\C* >= Vcc - Vdd + boVSc, (2.23) 

where 

VXA=v*il3)+v*{2A)-VAA 
V^B =v*{U) + v*{23)-Vbb 

VSc =v* (12) + v* (34) -Vcc. 

In the following fit bo is expected to be of the order unity and, indeed, fixing 6o = 1 is found to be a good assumption. 

The two-quark potentials v{ij) are taken to be more elaborate than the three term form of Eq. 2.8. They are fitted 
to the lattice data using [12] 

255 045 

v{rij) = 0.562 + 0.0696rjj- - —. (2.24) 

Similarly, the excitation with £J„ symmetry is fitted by 

/ -N 4.24 3.983 

Av = v*iij)-v{ij) = 2- + ^-- 2.25 

y ' ij ' ij 
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The extra terms containing r^^^ and r^^^ are purely for numerical reasons and ensure that the fitted values of v{ij) 
and v*{ij) are, on average, well within 1% of the lattice values for all > 2. 

e) There are two types of off-diagonal element 

< A*\V\B* > and < A\V\B* > . (2.26) 

Unlike Vab in Eq- 2.6, there is now no guide from the onc-gluon exchange limit of Eq. 2.5. Therefore, a further 

assumption is necessary for the form of these matrix elements. Several different forms were tried, but the most 
successful seems to be the one where the matrix elements needed in [V — {E + Vaa)N] = are expressed as follows 



<A*\V+^Vdd\B* >-- 



■y [Vabcd + coVXbc] =< B*\V + yVdd\A* > 



< A*\V+^Vdd\C* >= [Vcabd + coVSab] =< C*\V + ^Vdd\A* > 



< B*\V + ^Vdd\C* >=-^[Vbcad + cqV^ca\ =< C*\V + ^Vdd\B* >, 



(2.27) 



where 



Vabcd = Vaa + Vbb — Vcc — Vdd , Vcabd = Vcc + Vaa — Vbb — Vdd, 

Vbcad = Vbb + Vcc — Vaa — Vdd , 
Vabc = VXa + Vbb ~ ^cc> 

VSc + VXa - V^B , V^cA = V^B + VSc " V 



V 



CAB 



AA- 



These are seen to be simply the generalisation of Eq. 2.6 to the interaction involving excited gluon states. Also in 
analogy with bo in Eqs. 2.24, cq is a free parameter - hopefully of order unity. 

Likewise, 



<A\V+^Vdd\B* >=-y 



Vabcd + ao 



V 



ABC 



=< A*\V + ^Vdd\B > 



(2.28) 



where, again in analogy with bo in Eqs. 2.24, ao is a free parameter - hopefully of order unity. However, it should be 
added that this hope is even less well founded than the one for cq. 

In the special case of regular tetrahedra Vdd = Vaa = Vbb = Vcc and V reduces to the form 



(2.29) 



Vaa 

-!VaaI2 
-!VaaI2 


-/Vaa/2 

Vaa 
-!VaaI2 


-fVAA/2 
-fVAA/2 

Vaa 




-rva/2 
-rva/2 


-rva/2 



-rva/2 


-.rva/2 - 
-rva/2 






-rK/2 
-rK/2 


-rvai2 



-rvj2 


-rva/2 
-rvj2 




Vb 

-rK/2 
-rK/2 


-rvj2 

Vb 

-rvci2 


-rvj2 
-rv,i2 

Vb . 



where Va = Vaa + aoV^A/2 > H = Vaa + ^qV^a . K = Vaa + cqV^^- ^ '^it^ ^ geometries 



■ 1 

-//2 
-//2 


-//2 
1 

-//2 


-//2 
-//2 
1 




-r/2 
-r/2 


-r/2 



-r/2 


-r/21 
-r/2 






-r/2 
. -r/2 


-r/2 



-r/2 


-.r/2 
-r/2 




1 

-r/2 
-r/2 


-,/72 
1 

-r/2 


-r/2 
-r/2 
1 



(2.30) 
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The full 6x6 matrix [V — {E + Vaa)'^] now factorizes into three 2x2 matrices, two of which are identical - giving 
the observed degeneracy. These have the form 



[v-{E+VAAm 

whereas the third 2x2 matrix is 

[V-{E + VaaM = 



-E{l + f/2) -f<^(E-Va)/2 
-r{E - Va)/2 -E{1 + /72) + Vb + rVc/2 



0, 



-E{l-f) riE-Va) 

r{E - Va) -E{i - n + Vb- rvc 



= 0. 



(2.31) 



(2.32) 



In this case the problem reduces to solving two quadratic equations for E. However, away from the regular tetrahedron 
the complete 6x6 matrix needs to be treated directly. 



III. RESULTS 



This section is in two parts. In subsection III A the results for the 6 basis state model at /3 = 2.4 are given. There 
it is shown that only 5 of the possible 8 parameters have any significant influence. Also it is pointed out that the 
restriction to the 2 basis state version of subsection II A is distinctly inferior and that the further restriction to only 
two-body interactions (i.e. kf = in Eq. 2.11) is for most geometries very poor. In subsection IIIB, it is shown that 
the parameters of the model do not change significantly as /3 increases i.e. as the continuum is approached. Therefore, 
the parameters extracted at /3 = 2.4 could be used directly in, for example, a Resonating Group calculation of a 
four-quark model of meson- meson scattering as in Ref . [4] . 



A. The 6 basis state model at /3 = 2.4 



In Ref. [1] four quark energies have been extracted for a variety of geometries using a 16^ x 32 lattice with P = 2.4. 
Prom these energies, only one hundred - distributed over six geometries - are selected for fitting. Configurations 

containing flux links of less than two lattice units were not inchided, since they have strong lattice artefacts. In detail, 
we use 15 Tetrahedra (T), 6 Squares (S), 12 Rectangles (R), 4 Quadrilaterals (Q), 9 Non-Planar (NP) and 4 Linear 
(L). Only the lowest two energies (£^o,i) from the lattice simulation are used. But, as mentioned earlier before Eq. 2.8, 
the values of these two energies are more or less the same from the 2x2 and 3x3 simulations. In the latter case a 
third energy {E2) is in fact also available. However, it is not expected to be very reliable due to the higher excitations 
it contains. Its main purpose was to improve the estimate on Eqa by reducing the contamination from the higher 
states. The stability of (-Eo.i) is understandable, since they are the lowest states of given symmetry. For example, 
with squares their wavefunctions are simply i;{Eo) ^ (A + B)/V2 and ij{Ei) = {A- B)/V2. 

Before commencing a fit the size of the errors on the above data must be decided. The lattice simulations, through 
the boot-strap method, do indeed produce errors - statistical ones. However, some estimate must also be added for 
systematic errors. How this is done is somewhat subjective. Here the prescription is to assume all errors must be at 
least 0.005 and, also, at least 10% of the eigenvalue itself. The former corresponds to about 10% for the largest values 
of Eq and amounts to about 8 McV. 

The above 100 pieces of data are fitted with Minuit - the Migrad option. In the model there are, in principle, 
eight parameters: kf in Eq. 2.11, fi,f2,ka in Eq. 2.22, kc in Eq. 2.21, 60 in Eqs. 2.23 and ao,Co in Eqs. 2.28. The 
strategy for finding the optimal values of these parameters is as follows. Firstly, /f is fixed at zero and bg fixed 
at unity, as discussed after Eqs. 2.22 and 2.23. Furthermore, preliminary minimizations indicate that /Cc = is a 
good approximation and that ao, /2 are strongly correlated. This latter point is understandable, since for regular 
tetrahedra, where the 6 basis model is needed, the Vabcd, Vcabd etc. are zero, so that 

< A\V+^Vdd\B* >= -bJSSexp{-b,kaS)'^^^^SC 

i.e. only the combination ao/2 arises. To avoid this oq is, at first, fixed at four, the value that yields essentially the 
smallest x^/dof. This results in the parameter values shown in the first row of Table I. Releasing kc and bo to be free 
parameters gives the results in the second row. Finally, releasing also ao gives the third row. 
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Table I shows the following points: 

a) The value of kc is consistent with zero. This appears to be a common feature with all fits. This suggests that 
/° « 1, i.e. the excited configurations interact amongst themselves simply through 2-body potentials without the 
4-quark effect introduced when kc ^ 0. 

b) The second row shows that the expected value of 60 = 1-0 persists even when &o is allowed to vary. 

c) The third row shows that there is an instability between and ao, with both having errors that are Ri 25% of the 
parameter value. This instability also leads to a value of 60 that is much less than the expected value of unity. 

d) In the first two rows fco « S/f - in qualitative agreement with the approximate relationship /cq = 4/2 derived in 
the Appendix. 

e) The fourth row comes from Ref. [13], which contained a preliminary fit to the above data, but with a model in 
which the off-diagonal matrix elements of Eqs. 2.27, 2.28 had a slightly different form. The outcome, however, is 

qualitatively rather similar to the present results. 

In order to see from where the x'^ originates, Table II shows the contributions of the average 

N(G) 

= ]v(G) ^ - 

from each of the 12 types of data - i.e. from Eq,Ei for the six geometries {T,S,R,Q,NP,L). Here N{G) is the 
number of data points for geometry G, Mi are the model results and Ei ± AE^ are the lattice data. 

Another alternative for /" - not supported by the Appendix - is 

r = in + n^) eM-^sk'^s). (3.1) 

This gives a comparable x^/d.o.f. of 1.10, but has the inferior feature that /g is not consistent with zero and so in 
violation of the condition before Eqs. 2.22. 

For comparison, Model la with its one free parameter kf in Eq. 2.11 gives from Minuit a x^/d.o.f. of 3.16 for 
kf = 0.571(12). As seen from Table II most of this x^/d.o.f. comes from the Tetrahedral geometry - especially from 
the regular Tetrahedra. This was already anticipated at the end of Section II C. 

Another comparison is the extreme limit of / = 1, where only two-quark potentials are used in Model I. Here all 

the parameters are fixed, so that this is not a minimization. From Table II it is seen that this gives a very poor 
representation of the data - especially for Ex. As expected, the larger the size of the geometry the larger the 
contribution to x^/d.o.f., since those are the cases that need most the smaller values of /. For example, with a 
6x6 square the model gives i?o,i = —0.124,0.3721, whereas from the lattice simulation the corresponding energies 
are —0.028(5), 0.065(7). This comparison is the basis of our statement in the introduction that multi-quark models 
containing only the standard two-quark potential of Eq. 2.8 or even the more elaborate form in Eq. 2.24 seem to have 
no justification for the cases considered here. 



B. The continuum limit 

The above results are for (3 = 2.4, which corresponds to a lattice spacing of a ~ 0.12 fm. But, in practice, if the model 
is to be used, for example, in some type of resonating group approach [4], then it is the continuum limit that is needed. 
Such an extrapolation would, in principle, require repeating the above lattice calculations with all 6 geometries for a 
series of increasing values of /?. This we do not plan to do here. Instead we exploit the conclusions of Ref. [11]. There 
it was shown that the four-quark binding energies already scale at /3 = 2.4, evcnthough the total two- and four-quark 
energies do not scale separately. Therefore, here we simply assume that, if some binding energy E{(}) would be scaled 
to the /3 = 2.4 lattice by the transformation pE[f3,r{l3)/ p] E{2A), where p = a(2.4)/a(/3), then we would find 
that this scaled £'(2.4) equals the binding energy E{r) we have calculated directly on the f3 = 2.4 lattice. Here the 
r{(3) and r are distances in the lattice units of the two different lattices. However, to completely specify the model 
at some other value of /? also needs a knowledge of v{ij,(3) and Av{ij,f3) corresponding to the f3 = 2.4 expressions 
in Eqs. 2.24, 2.25. In Ref. [11] it was shown that v{ij,(3) did not scale. But this is not a problem, since v{ij,f3) is 
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documented in the literature for other values of /? - see Refs. [14-16]. These can then be scaled in the same way 
as the E{f3) to give v{ij,2A). However, as found in Ref. [11] we do not expect these 2.4) to equal the v{ij) in 
Eq. 2.24. It remains to specify Av{ij,(3). Here we simply assume that it scales, so that the values at /? = 2.4 can 
be used directly in precisely the same way as the binding energies E{(3). The rational behind this last assumption is 
twofold: 

i) In the model it is the differences Av and v{ij,(3) — v{kl,(3) that enter as in Eqs. 2.23, 2.27 and 2.28. Therefore, 
some of the lack of scaling in v* {ij, (3) and v{ij, (3) may well cancel - as in the four-quark binding energies. 

ii) Any lack of scaling that remains will be reflected in a lack of scaling of the ao, 5o, co multiplying the Af and y*'s. 

The actual values of (3 used are 2.5 [14], 2.5 (denoted by 2.5'), 2.635 and 2.74 [15] and 2.85 [16]. In the last case, 
the authors present two potentials[y(i?, 5) and V^(oo)] that are appropriate for interquark distances R greater than 
6a(2.85) « 0.17 fm. Their own preference is V{R, 5) - the one we use. 

The results arc shown in Tabic III. The main observation is that, of the free parameters, only fc/ appears to be a 
significant function of /3 - but even there all values are within the error bars. Therefore, the final version of the model 
has only four independent parameters. 

A similar continuum extrapolation was carried out in Ref. [11]. However there are several major differences compared 

with the above analysis: 

a) Much less data for each value of /? was included in the fits, since only Squares and Tilted Rectangles were analysed. 
Furthermore, only the Eqs were fitted. This means that, for each /?, there were only 17 pieces of data compared with 
the present case of 100 pieces covering 6 geometries with both Eq and E\ . 

b) All the two- and four-quark energies were simulated at five values of (3 and interpolated to the same physical sizes. 
As scaling was found for the binding energies, they are here obtained directly from the (3 — 2.4 values. 

c) In order to study scaling and the effect of lattice artefacts, the errors on the data were taken to be purely statistical 
and so were much smaller than the present ones. 

d) The model analysis was carried out with the 2x2 model. Since this results in a quadratic equation for the energies, 
the step Eq f could be done analytically. After this the parameters kf, kp, /o in Eq. 2.13 were extracted by a 
fitting procedure. 



C. A comparison with other work 

A fit with a potential model to some of the above four-quark binding energies was attempted earlier in Ref. [17]. 

However, there are several important differences compared with the present fit: 

1) The model contained only 2 x 2 or 3 x 3 bases that were essentially constructed from A,B,C. 

2) The Tetrahedral data was not included. 

3) Most of the data was for P = 2.4. Only squares included that from f3 = 2.5. 

4) The overlap factor is that in Eq. 2.13 but with /g = 1. 

For the four geometries S,R,Q,L the area 5 is the same as that used here in Eq. 2.12 and for geometry NP a curved 

surface is used. However, the perimeters P are determined by the underlying lattice structure of the geometry and 
not the usually expected perimeter. For example, for Quadrilaterals, P is the perimeter of the rectangle enclosing 
the quadrilateral. The reason for this choice is an attempt to take into account lattice artefacts, which are present 
for the smallest configurations. These artefacts are more important in Ref. [17], since configurations containing fiux- 
lines of only one lattice link are included. In the present work the cut-off is at two lattice links. For Squares the 
outcome in [17] was kf = 0.296(11) and kp = 0.080(2), compared with kj = 0.571(12) using Model la above in 
Subsection II A. The difference between the values of these area constants is because the perimeter term, which is 
introduced to make artefact corrections for small configurations, is still important for large configurations, where 
rotational symmetry is restored and their presence not needed. For example, with a square of side d — 5, the two 
terms in the exponent of Eq. 2.13 are comparable, even though the lattice data for this square is -0.042(5) and 
-0.047(5) for the corresponding tilted square. In the lattice simulation the basis states have a similar form compared 
with the ones used in [17] but with the major difference that they are fuzzed to such an extent that these underlying 
basis states get transformed. Therefore, in [17], for larger configurations it would be necessary to include many other 
terms with different perimeters, so that on the average the usual perimeter for a particular geometry would probably 
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be more appropriate. This latter choice was tested in Ref. [11], where it was found to reduce in the continuum hmit 
to the form of Model la. Thus the perimeter term seems to measure lattice artefacts also when the physical perimeter 
is used. 

In Ref. [10], the model in Ref. [17] is extended to include the tetrahedral data by using an 8 x 8 basis, where now all 
the states are constructed in terms of the underlying lattice. This, therefore, suffers from the above problem that the 
results are always dependent on the lattice even for large configurations, where rotational invariance is restored. 



IV. CONCLUSIONS 

In this article an interquark potential model is constructed than can explain, on the average, the energies of a series 

of four-quark systems, in which the four static quarks are in any one of six geometries. The input energies used as 
data for this construction [1] were calculated on a 16'^ x 32 lattice for four static SU(2) quarks with f3 = 2.4. Two or 
three basis states were used in the lattice simulation enabling reliable estimates to be made of the lowest two energy 
eigenvalues (-Eo,i)- In aUj 100 pieces of data were considered suitable for confronting the model. Other pieces were 
rejected, if they involved configurations suspected of having strong lattice artefacts (e.g. with states constructed from 
less than two lattice links) or had vanishingly small binding energies (e.g. elongated rectangles with r ^ d). 

The full model utilized 6 basis states A, B, C, A* , B* , C* and in its most general form had eight parameters. However, 
in practice, only 4 of these (kf, ka, f2 and cq) need be considered as free when fitting the data. 

The parameters that are, perhaps, of most interest are those connected with the ranges of the various interactions, 
namely, kf,ka,kc- Here we define 'range' as rf^a,c = y^jbjzj^c- In Model la, where ka.c arc effectively infinite, 
we get fc/(/a)=0.57(l) i.e. rf{Ia) = 5.0 in lattice units. However, in Model III as the excited states A*,B*,C* are 
introduced, the interaction between the basic states A, B, C decreases by raising kj to 1.51 giving rf{III) = 3.1. But 
at the same time this loss of binding by the direct interaction between A, B,C' is compensated by their coupling to the 
A*,B* ,C* states. This coupling is found to have about the same range = 5.1 as rf{Ia) above, whereas the direct 
interaction between the A*,B*,C* states seems, in all fits, to be satisfied with simply a two-quark description without 
any four-quark correction (i.e. kc=0). The observation that rf{Ia) w suggests that the energy density has a range 
dictated by the longest range available - namely ra- Therefore, when the A* ,B* ,C* states are not explicitly present, 
as in Model la, the only available range rf{Ia) has to simulate the role of ra- In the binding energies the contributions 
from the A*,B*,C* states rapidly dominate over those from the A,B,C states. For example, with squares of side 
R, the A, B, C states contribute only 85, 40, 10% to the binding energy for i?=2,4,6 respectively. Of course, at the 
largest distances (« 0.7 fm) the quenched approximation is expected to break down and the role of quark-pair creation 
to become important. We, therefore, come to the following scenario for the four-quark interaction. At the shortest 
distances, upto about 0.2 fm, perturbation theory is reasonable with the binding being given mainly by the A, B, C 
states interacting simply through the two-quark potentials with little effect from four-quark potentials. However, for 
intermediate ranges, from about 0.2 to 0.5 fm, the four-quark potentials act in such a way as to reduce the effect of 
the A, B, C states so that the binding is dominated by the A*, B*,C* states, which now interact amongst themselves 
again simply through the two-quark potentials with little effect from four-quark potentials. At larger ranges quark-pair 
creation can no longer be neglected. 

The next step in this work would be to carry out a Resonating Group calculation in the spirit of Ref. [4]. This 
would then test the universality of the above model by extending it from the 6 geometries discussed here to general 
four-quark geometries. 
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The justification of Eq. 2.22 



In the case of regular tetrahedra a simple model for the ^4(13, 24), A*{13, 24) states is as follows: 

A{13,24) = N^{ab + cd)^{ef + gh), 



5(14,23) = N-^{aj + ih)-^{ek + Id), 



A* {13, 24) = N*^ {ah - cd) (e/ - gh) , 



B*(14, 23) = N*-^{aj - ih)-^{ek - Id), (5.1) 
V 2 v 2 

where a, b, c.k, I are links along the sides of the cube defining the regular tetrahedron and A'^ and A''* are normalisation 
constants. The loops a, b, c, d and e, /, g, h are taken to be in the xz plane and the loops a,j, i, h and e, k, I, d are in 
the yz plane. 

The above definitions guarantee < A(13, 24)^* (13, 24) >= etc. and give for the two normalisations: 

< AA >= N'^[l + I + S{xz) + S{xz)f/A i.e. N = 1/[1 + S{xz)] and 

< A* A* >= 7V*2jl + 1 - S{xz) - S{xz)]'^ i.e. N* = 1/[1 - S{xz)], 

where S{xz) is a Rx R Square loop of links in the xz plane along the sides of the cube i.e. perimeter is of length AR. 
Other loops needed in this discussion are as follows: 

The loop L is two connected faces of the cube in an L-shape i.e. perimeter is of length 6R, 
C is three connected faces of the cube giving a corner i.e. perimeter is also of length 6R 
T is three faces of the cube giving a "table" i.e. perimeter is of length 8R 
With this notation the following overlaps can be written down: 

< A(13, 24)5(14, 23) >= -^ji^S + 2T + AC + 8L]/(1 + S)^ = -f/2, 
where the initial factor of -1/2 is the usual colour factor giving 

^ _ [25 + 2T + 4C + 8L] 



4(1 + 5)2 



Likewise, 



[2S + 2T + 4C-8L] [2S + 2T-4C] 

~ 4(1-5)2 ^'"'^ ^ ~ 4(1-52) • 

A rough estimate of the relative importance of S,L,C,T can be seen as follows: In a Wilson loop, the further two 
combinations of Unks (e.g. LI, L2) are separated in Euclidean time, the less is their correlation < L1.L2 >= exp{—VT). 
Similarly, spatial links are less correlated the further they are apart. 

For the Square (5) let side 1 have a factor 1 and the two adjacent sides factors of p and the fourth side p^. This 
makes 5 = DA.p.p.p^ = p*, where D is some constant and p some correlation factor to adjacent sides with p <1 
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For the L-shape let the correlation factors have an extra q, when changing direction out of the original side of the 
cube and an extra r for the final side. In all L = D.l.p.pq.pqr.pq.p = Dp^q^r when moving around the perimeter of 
the Z/-shape. 

For the C-shape the same factor arises as in the i-shape. This suggests L = C. 

For the T-shape the last set of links is close to original ones, so that the factor on the last link is as in the S case. 
In all T = D.l.p.pq.pqs.p^.pqs.pq.p = Dp^q^s^. 

For large R we expect p,q,r,s <^ 1 which suggests S ^ {L or C) ^ T. 
For small R we expect p,q,r,s ^ 1 i.e. S" (L or C) « T w 1 
One way to parametrize these two limits is to take 

S = exp{—WR'^) pa (1 — WR^) since it involves a single square of area R"^. 
L = ex.p{—W2R'^) w (1 — 2WR'^) since it involves two squares of area R"^. 
T = exp(— W^3i?^) « (1 — SWR^) since it involves 3 such squares. 

Keeping only the terms upto 0{B?) in S, L, C, T gives simply / = 1 — WR^, = /° = and is model independent. 
Whereas the addition of the terms from expanding further the exponentials gives f" = WR'^/A and p = 1/2. 

This last result clearly shows that the 0{R'^) terms are not so simple, since we expect /'^ ^ 1 as i? ^ - as in the 
case of /. Similarly these arguments would suggest kf = 4/2 . However, the R^ dependence of /" is probably correct 
and so for general geometries we assume f" oc Area as Area — > 0. For the case of regular tetrahedra, the appropriate 
area from Eq. 2.12 is VSR'^. 
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kf 




/I 


he 


ao 


bo 


Co 


xVdof 


Fix 4 


1.51(8) 


0.55(3) 


0.51(2) 


ao 




LO 


3.9(3) 


1.03 


Fix 2 


1.50(8) 


0.57(3) 


0.54(7) 


0.02(5) 


4.0 


1.06(.32) 


4.3(1.3) 


1.05 


Fix 1 


1.61(11) 


0.59(5) 


0.28(6) 


0.00(11) 


:5.4(8) 


0.19(2) 


0.78(6) 


0.90 


Ref. [Vr 


1.2.-)(()) 


()..") 1(11) 


().()()( 1) 


0.01(20) 


1. IC!) 


2.2(()) 


8.0(1) 


1.08 



TABLE I. The values of the parameters defining the interaction and fixing = 0. Those 
numbers underlined are kept fixed. 



Model 


Model III 


Model la 


Model I(f= 


=1) 


Geometry 


Xa{Eo) 


Xa{Ei) 


Xa{Eo) 


Xa{Ei) 


XAiEo) 


XAiEi) 


T 


0.9 


0.7 


8.5 


7.6 


15.8 


141.5 


S 


0.2 


1.9 


1.2 


1.5 


102.4 


659.6 


R 


1.1 


0.9 


1.1 


0.5 


79.2 


527.2 


Q 


2.7 


1.2 


1.0 


0.3 


0.8 


6.2 


NP 


1.0 


0.2 


1.7 


0.2 


0.6 


16.9 


L 


0.3 


2.6 


0.3 


2.6 


0.3 


2.6 



TABLE II. The contributions Xa(^o,i) from the two states {Eq, Ei) in each of the six geometries 
(T, S, R, Q, NP, L) to the total x^/d.o.f. of 1.03 for Model III, 3.16 for Model la and 144 for Model 
I(/ = l)- 



Paraiiii'tiT 


J = 2.1 


J = 2..") 


J = 2..r 


:S = 2.():!5 


J = 2.71 


:S = 2.85 


kf 


1.51(8) 


1.48(8) 


1.50(8) 


1.49(8) 


1.45(8) 


1.43(8) 


ka 


0.55(3) 


0.56(3) 


0.56(3) 


0.56(3) 


0.57(3) 


0.57(3) 


/2 


0.51(2) 


0.51(2) 


0.51(2) 


0.51(2) 


0.51(2) 


0.51(2) 


Co 


3.9(3) 


3.9(3) 


3.9(3) 


3.9(3) 


3.9(3) 


3.9(3) 


x7d-o.f. 


1.03 


1.03 


1.11 


1.09 


1.15 


1.07 



TABLE III. The values of the parameters defining the interaction fixing f^ = kc = 0, bo = 1.0 
and ao = 4.0 for a series of /3's. 
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